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We perform a timc-dcpondont and time-integrated angular analysis of the decays 
(pK* {892f, ipK^iUSOf, and (p{KTr)o° with the final sample of about 465 million SB pairs recorded 
with the BAR4R detector. Twenty-four parameters are investigated, including the branching frac- 
tions, CP-violation parameters, and parameters sensitive to final-state interactions. We use the 
dependence on the Kiv invariant mass of the interference between the scalar and vector or tensor 
components to resolve discrete ambiguities of the strong and weak phases. We use the time-evolution 
of the B ipK%iT^ channel to extract the CP-violation phase difi^erence A(/)oo = 0.28 ± 0.42 ± 0.04 
between the B and B decay amplitudes. When the B — » ipK^'K^ channel is included, the fractions 
of longitudinal polarization /z, of the vector-vector and vector-tensor decay modes are measured to 
be 0.494 ± 0.034 ± 0.013 and 0.90lto:o58 ± 0-037, respectively This polarization pattern requires 
the presence of a positive-helicity amplitude in the vector-vector decay from a currently unknown 
source. 
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PACS numbers: 13.25.Hw, 13.88.+e, ll.30.Er, 12.15.Hh 



I. INTRODUCTION 

Charge-Parity (CP) symmetry violation has been rec- 
ognized as one of the fundamental requirements for pro- 
ducing a matter-dominated universe [Ij and therefore it 
has played an important role in understanding funda- 
mental physics since its initial discovery in the K meson 
system in 1964 0. A significant CP- violating asymme- 
try has been observed in decays of neutral B mesons to 
final states containing charmonium, due to interference 
between B^-B^ mixing and direct decay amplitudes 
It has been established 4] that the CP-violating decays 
of the meson are due to CP violation in decay ampli- 
tudes, as well as in K^-K'^ mixing, and this kind of "di- 
rect" CP asymmetry in B decays has also been observed 
recently The CP asymmetries are generally much 
larger in B decays than in K decays Q because they di- 
rectly probe the least flat Unitarity Triangle constructed 
from the Cabibbo-Kobayashi-Maskawa (CKM) matrix el- 
ements Vij , which relate weak and flavor quark eigen- 
states [tI . This Triangle reflects the unitarity of the CKM 
matrix, and two of its angles are the phase differences of 
its sides on the complex plane: a = arg(— Vj^jV^j/T^dKTb) 
and /3 ^ arg(-K,K;/y,,y;j. Due to the large CP- 
violating effects, B decays provide an excellent testing 
ground of fundamental interactions. 



B 




d 



K 

d 



FIG. 1: Penguin diagram describing the decay B ipK^*\ 
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The CP-violating effects observed to date are self- 
consistent within the Standard Model with a single com- 
plex phase in the CKM mechanism However, this 
mechanism alone is believed to be insufficient to produce 
the present matter-dominated universe. Therefore, it is 
important to search for new sources of CP-violating in- 
teractions. While direct access to new fundamental par- 
ticles and interactions may be beyond the energy reach of 
operating accelerators, one can look for them in virtual 
transitions. New particles in virtual transitions, includ- 
ing but not limited to supersymmetric particles Q , would 
provide additional amplitudes with different phases. De- 
pending on the model parameters, sizable CP-violating 
effects, either B^-B^ mixing-induced or "direct" , could 
be observed in pure penguin modes which involve virtual 
loops as in the example shown in Fig. [T] Some of the 

first observed gluonic penguin decays, B rj' K P 

and B (pK^*^ [3, remain promising channels in 
which to look for new physics. The latter type of decay 
is illustrated in Fig. [T] and is the focus of this paper. For 
example, comparison of the value of sin 2/3 obtained from 
these modes with that from charmonium modes such as 
[Tl| , or measurement of direct CP vi- 
olation, can probe new physics participating in penguin 
loops [12]. 

The (V-A) nature of the weak interaction leads to left- 
handed fermion couplings in interactions with W bosons, 
such as those shown in Fig. [TJ Combined with helicity 
conservation in strong interactions and spin-flip suppres- 
sion of relativistic decay products, this leads to certain 
expectations of the spin alignment in weak B meson de- 
cays to light particles with spin, such as B ip K * p^ . 
However, the large fraction of transverse polarization in 

the 'i?' ^ 'jc'*(892) decay measured by BaBAR fTl| and 



by Belle 15| indicates a significant departure from the 
naive expectation of predominant longitudinal polariza- 
tion. This suggests the presence of other contributions to 
the decay amplitude, previously neg lected, either within 
or beyond the Standard Model [16]. The presence of a 
substantial transverse amplitude also allows the study of 

CP violation in the angular distribution of B <p K * 
decays, an approach complementary to either mixing- 
induced or yield asymmetry studies. Polarization mea- 
surements in B decays are discussed in a recent re- 
view [13, di- In Table U we list BaBAR's recent mea- 
surements of the branching fraction and longitudinal po- 
larization in the decays 's' ^ (p ^K^[j*^ [H, [13, HU, IH, HI . 

Measurements in B p K * decays have also revealed 
a large fraction of transverse polarization ^J] . 

In this analysis, we use the final sample of about 465 
million T(45') BB pairs recorded with the BaBAR de- 
tector at the PEP-II asymmetric-energy e+e^ storage 
rings at SLAC. We employ all these techniques for CP- 
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TABLE I: BABAR's recent measurements of the branching 
fraction B and longitudinal polarization fraction in the de- 
cays B — > (pKj*\ The spin J and parity P quantum numbers 
of the /Cj*' mesons are quoted. The upper limits are shown 
at the 90% confidence level. For a complete list of all observ- 
ables in each analysis see the Refs. listed. Results indicated 
with f are superseded by this analysis. 



Mode 


J 


Ref. 


B (10 ) 


/i 




0- 


fl9] 


8.4lij ±0.5 


1 


ipK+ 


Q- 


fl9] 


lO.Olj;:^ ± 0.5 


1 


<y3A'o*(1430)'' 


0+ 


[20jt 


4.6 ±0.7 ±0.6 


1 


(pA'o*(1430)+ 


0+ 


[2^ 


7.0 ±1.3 ±0.9 


1 


95 a:* (892)° 


1" 


[20jt 


9.2 ±0.7 ±0.6 


0.51 ± 0.04 ± 0.02 


ipK*{m2)+ 


1" 


[22] 


11.2 ±1.0 ±0.9 


0.49 ±0.05 ±0.03 


v5is:*(i4io)+ 


1" 


[2^ 


< 4.8 






1~ 


[23] 


< 3.5 




(^^-1(1270)+ 


1+ 




6.1± 1.6± 1.1 


r, AC+0.12 +0.03 
"•^"-0.13 -0.07 


if>Ki{im)+ 


1+ 


[21] 


< 3.2 




(^7^1(1430)° 


2+ 


[20jt 


7.8 ± 1.1 ±0.6 


0.85lo:o? ± 0.04 


(^7^2* (1430)+ 


2+ 


[2^ 


8.4 ± 1.8 ±0.9 


0.80lHo ± 0.03 


^3^-2(1770)+ 


2" 


[2^ 


< 16.0 




(^7^2(1820)+ 


2~ 




< 23.4 






3" 


[23] 


< 2.7 






4+ 


[23j 


< 15.3 





violation and polarization measurements in the study 

of a single B-decay topology 5° — !■ (p{Ktt). Over- 
all, 27 independent parameters sensitive to CP violation, 
spin alignment, or strong- or weak-interaction phases 
describe three decay channels (twelve in either vector- 
vector or vector-tensor and three in vector-scalar de- 
cays), which leaves only one overall phase unmeasur- 
able. The three channels in our amplitude analysis are 

W° ^K*(892)0, V? '^'2(1430)°, and ipCK'K)^. The 

latter contribution includes the K q(1430)° resonance to- 
gether with a nonresonant component, as measured by 
the LASS experiment [25|. While we describe the analy- 
sis of these three neutral- i? meson decays, this technique, 
with the exception of time-dependent measurements, has 
also beeii applied recently to the charged- i? meson de- 
cays [H . 

We use the time-evolution of the _B " ^ ipKgi:^ chan- 
nel to extract the mixing-induced CP- violating phase dif- 
ference between the B and B decay amplitudes, which is 
equivalent to a measurement of sin 2/3 to a good approxi- 
mation. With the B ° ^ (fK^n^ channel included, the 
fractions of longitudinal and parity-odd transverse ampli- 
tudes in the vector- vector and vector-tensor decay modes 
are measured. We use the dependence on the Kt: invari- 
ant mass of the interference between the scalar and vec- 




mK.(GeV) 



FIG. 2: Invariant Kn mass distribution from the B 
(fiK'^TT^ analysis from Refs. [13, HI]. The solid (dashed) 
line is a projection of the signal-plus-background (background 
only) fit result. The narrow charm background peak at 1.865 
GeV comes from decays to Ktt and is not associated with 
ipK^Tv^ production. The arrow indicates the mass range con- 
sidered in this analysis. 



tor, or scalar and tensor components to resolve discrete 
ambiguities of the strong and weak phases. Using either 
interference between different channels or B^-B'^ mixing, 
we measure essentially all 27 independent parameters ex- 
cept for three quantities that characterize the parity-odd 
transverse amplitude in the vector-tensor decay, which is 
found to be consistent with zero. 



II. ANALYSIS STRATEGY 

Earlier studies of -B ° — > (pK^ir^ decays by the BaBAR 
collaboration [l^, H^l indicate the presence of three sig- 
nificant Kit partial waves: (Kir)'^'^ (spin J = 0, includ- 
ing the resonance i4:o*(1430)''), K*(892)0 (J = 1), and 
i4:|(1430)° (J = 2). These correspond to the following 
decays, with the number of independent amplitudes char- 
acterizing different spin projections given in parenthe- 
ses: ^ <^('^'7r)*o (one), *i?'*(892)0 (three), and 

(y9 if 2(1430)° (three). No significant contribution from 
other final states has been found with Ktt invariant mass 
rriKTv up to 2.15 GeV [HI, HI]. See Fig.[2]for an illustration 

of the B ^ ^ ipK^TT^ contributions. Therefore, we limit 
our analysis to the mass range niKTr < 1.55 GeV with- 
out any significant loss of i? ^ ipK^-K^ signal through 
charmless Kt: resonant or nonresonant production. 

There has been no extensive study of the B ^ ^ 

ipKgTT^ decay, except for the study of B ^ 

ipK* (892)'^ [l3|. However, due to isospin symmetry 
of the if^Tr" and i^^'^Tr^ systems, the same amplitude 
composition is expected in the ipK^n^ and ipKgiT^ final 
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states. We do not expect any charmless resonance struc- 
ture in the tpK^ or i/jtt^ combinations, while we veto the 
charm resonance states, such as D^-j (pir^. 

It is instructive to do a simple counting of the ampli- 
tude parameters in B ^ ^ (p K tt decays with the three 
Kn spin contributions discussed above. With seven in- 
dependent Ajx complex amplitudes for B decays and 
seven Aj\ amplitudes for B decays, we could construct 
28 independent real parameters. Here J refers to the 
spin of the Kn system and A to the spin projection of 
the if meson onto the direction opposite to the B me- 
son flight direction in the ip rest frame. However, one 
overall phase is not measurable and we are left with 27 
real measurable parameters. Among these parameters, 

26 parameters have been or can be measured in the de- 
cay B (pK^TT^ [23|. Those are branching fractions, 
polarization parameters, strong phases, and CP asymme- 
tries. Some of the phases are extracted from the inter- 
ference effects between different modes. However, due to 
limited statistics some of the CP asymmetries were not 
measured in earlier analyses and we now extend those 
measurements. 

Finally, one parameter, which relates the phases of 
the B and B decay amplitudes, can be measured us- 
ing only the interference between decays with and with- 
out B'^-B^ mixing, to final states which can be decom- 
posed as CP eigenstates, such as pK^n^ . In Table [TTl all 

27 real parameters measurable with B ^ p K tt de- 
cays are summarized. These parameters are expressed in 
terms of the Aj\ and Aj\ amplitudes for pK^ir^ 
or pK^TT^ and B° pK~iT^ or pK^n^ decays. We 
also refer to a transformed set of amplitudes Ajq and 
Aj±i = ± Aj±)/^/2. The parameters in Table |TT] 
are expressed as six CP-averaged and six CP-violating 
parameters for the vector-vector and vector-tensor de- 
cays. The TT in the definitions of 4>±j and accounts 
for the sign flip A±j = —A±j if CP is conserved. The 
parameterization in Table [H] is motivated by the negligi- 
ble CP violation expected in these decays. Therefore, the 
polarization parameters specific to either B (superscript 

") or B (superscript "+") are the CP-averaged pa- 
rameters with small CP-violating corrections which are 
either multiplicative (for rates) or additive (for phases): 





= Bj-{l±Acpj)/2 


(1) 




= hj-il±A%j) 


(2) 


fh 


= Uj ■ (1 ± Ahpj) 


(3) 






(4) 




= (^xj ± A0ij + 1 ± 1 


(5) 




= ^oj±A(5oj. 


(6) 



distributions critical to separating different partial waves, 
then introduce interference effects between amplitudes 
from different decays, and finally discuss time-dependent 
distributions. 



A. Angular distributions 

We discuss here the angular distribution of the decay 
products in the chain B — > pK* (K'^ K~){Ktt) in- 
tegrated over time. First we look at the decay of a P 
meson only and leave the B for later discussion, which 
involves CP violation. Angular momentum conservation 
in the decay of a spinless B meson leads to three possi- 
ble spin projections of the p meson onto its direction of 
flight, each corresponding to a complex amplitude Ajx 
with A = or ±1. The three A values are allowed with 
the K* spin states J > 1, but only A = contributes 
with a spin-zero K*. The angular distributions can be 
expressed as functions of Hi — cosOi and <I>. Here 6i 
is the angle between the direction of the K meson from 
the K* Kit {Oi) or p KK [62) and the direction 
opposite to the B in the K* or </? rest frame, and $ is 
the angle between the decay planes of the two systems, 
as shown in Fig. [31 The differential decay width is 



dHi dH2 d^ 



AjxYf{ni,<^)Yf\-n2,0) 



(7) 



where Yj are the spherical harmonic functions, with J = 
2 for i^*(1430), J = 1 for K*{892), and J = for {Kn)*, 
including ArQ(1430). We do not consider higher values of 
J because no significant contribution from those states is 
expected. Only resonances with spin-parity combination 
P — (—1)'' are possible in the decay K* Ktt due to 
parity conservation in these strong-interaction decays. 

If we ignore interference between modes with different 
spins J of the Kn system in Eq. ([7|), then for each de- 
cay mode we have three complex amplitudes Ajx which 




In this Section we discuss further the method for the 
measurement of the relative phase, along with all the 
other parameters. First we review the angular distribu- 
tions, follow with a discussion of the Kt: invariant mass 



FIG. 3: Definition of decay angles given in the rest frames of 
the decaying parents. 
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TABLE 11: Definitions of 27 real parameters measurable with the B (pKn decays. Three resonance final states with spin 
J = 0, 1, 2 are considered in the Kn spectrum. The branching fraction B is calculated as a ratio of the average partial decay 
widths for (F) and B'^ (F) and the total width Ftotai where we neglect any difference in the and B'^ widths. This 
definition allows for differences between the _B° and _B° decay amplitudes, Aja and Ajx, as discussed in the text. 



parameter definition ipK^{U30) if K* (892) ifiK^{U30) 

Bj i(rj + rj)/rtotai 

1 (•] T ... |2 /v^l TT _ |2 I I /I _ |2 /.pl < _ |2\ 



fLj i(lA/o|7s|A/Ar + |A/o|7S|A. 



^CPJ 



f±j i(|Ajx|7s|A/A|' + |A/±|Vs|A,7a|') 

hj 5(ai-g(Aj||/A/o) + arg(ylj||/Ajo)) 

0XJ i(arg(Ijx/ljo) + arg(AjxM.jo) - vr) 

Soj i(arg(Aoo/^jo) + arg(^ooM.Jo)) 

Acpj (f J - rj)/(f J + Fj) 



A^PJ (|^.7x|7S|^JAr - |^Jx|7S|A/Ar)/(|^.7x|7S|^JAr + |A/x|7S|A 

A0IIJ i(arg(Aj||/A7o) - arg(4,j||/Ajo)) 

A0XJ |(arg(Ajx/Ajo) - arg(AjxMjo) - vr) 

A5oj i(arg(Zoo/Ajo) - arg(^ooMjo)) 

A0OO iarg(A)oMoo) 



JA 



J = 


J = 1 


J = 2 


Co 


Bi 


B2 


1 


hi 


/l2 


none 


f±i 


f±2 


none 


hi 


<t>\\2 


none 


0X1 


0X2 





<5oi 


5o2 


Acpo 


Acpi 


AcP2 





A" 

■^CPl 




none 


Acpi 


AcP2 


none 


A0II1 


A0II2 


none 


A0X1 


A0X2 





AJoi 


A5o2 


A 000 


none 


none 



appear in the angular distribution. We discuss interfer- where the functions fij (Tii, 7^2, are given in Ta- 
ence between different modes later in this Section. The ble IIIIl The a~j parameters are defined as: 
differential decay rate for each decay mode involves six 
real quantities a~j, including terms that account for in- 
terference between amplitudes of common J. 

J2o.-jUj{ni,n2,'^>), (8) 



"5.7 



=/l7 (9) 



|Aj|| p + |Aj^p _ |A,;+ip + |Aj_ip 

^|A,7aP 

|A,,||p-|A,,^|2 _ 5Re(A,,+iA}_i) 



= S|A,.P = ' I]|A,.p" = - ^^^^ 



" ^A^TP ^ ■ E|A,,p" " " ^^'^ ^ ■ ^^^^ ^^^^ 
Sm(A.,^A}||) Qm(A,,+iA}_,) / „ , — , , 

3^e(Aj||A}o) jRe(A,7+iA}o + A,7_iA}o) A—— — 

5r.(A,,A}o) _ ^^m(A,,, A}„ - A,_, A}o) ^ 7" . 



S|AjaP V2-S|AjaP 
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The above terms are specific to the decays and 
are denoted with the superscript "— ", as introduced in 
Eqs. HHSI). The angular distributions for the B° decays 



are described by the same Eq 
by a+ 

ing Ahy A and superscript 



but with a^j replaced 



jj, and with definitions given by Eqs. 



MT4j) . replac- 



by 



B. Mass distributions 

The differential decay width given in Eq. ([7]) is param- 
eterized as a function of helicity angles. However, it also 
depends on the invariant mass tokit of the Kn resonance, 
and the amplitudes should be considered as functions of 
'n^K-K- Without considering interference between differ- 
ent modes, as shown in Eq. ([8]), this mass dependence 
decouples from the angular dependence. Nonetheless, 
this dependence is important for separating different Kn 
states. The interference effects will be considered in the 
next subsection. A relativistic spin-J Breit-Wigner (B- 
W) complex amplitude Rj can be used to parameterize 
the resonance masses with J = 1 and J = 2 : 



Rj{m) = 



mjTj{m) 



(mj — m?) — imjT j{fn) 
where we use the following convention: 

cot 8 J 



sin (5,76 



id, J 



mjTj{m) ' 
The mass-dependent widths are given by: 

3 



TiM = El 



TOi 1 



m 1 + r^g^ \(7i 



TO2 9 



4 4 



(15) 



(16) 



(17) 



(18) 



where Fj is the resonance width, mj is the resonance 
mass, q is the momentum of a daughter particle in the 
resonance system after its two-body decay {qj is evalu- 
ated at m = uij), and r is the interaction radius. 

The parameterization of the scalar (_fir7r)g'^ mass distri- 
bution requires more attention. Studies of Kir scattering 
were performed by the LASS experiment [25| . It was 
found that the scattering is elastic up to about 1.5 GeV 
and can be parameterized with the amplitude: 



iSo 



where 



i?o(™) = sin(5oe 



So^ AR + AB, 



(19) 



(20) 



AR represents a resonant ifQ(1430)° contribution and 
AB represents a nonresonant contribution. The mass 
dependence of AB is described by means of an effective 
range parameterization of the usual type: 



cot AB 



aq 



(21) 



where a is the scattering length and b is the effective 
range. The mass dependence of AR is described by 
means of a B-W parameterization of a form similar to 
Eq. HID: 



cot AR = 



(22) 



TOoro(m) ' 

where ttiq is the resonance mass, and ro(TO) is defined as 



ro(r7i) = To 



mo / q 



m 
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(23) 



The invariant amplitude Mj{m) is proportional to 
Rj[m): 



Mj{m) cx — Rj{m) 

and can be expressed, for example for J = 0, as 
Mo(m) 



(24) 



„2jAB 



q cot AB — iq 
^oml/qo 



(toq — m^) — imoTo^m) 



(25) 



The resulting {KTr)Q'^ invariant mass distribution is 
shown in Fig. [4l along with the phase and distributions 
for the other resonances. The mass parameters describ- 
ing the three spin states in the m distribution are shown 
in Table HVl Measurements of the LASS experiment are 
used for the parameters of the J — contribution and 
for the interaction radius [1^, H^. The values of mo, Fq, 
a, and b used in this analysis are different from those 
quoted in Ref. 1251 due to better handling of the fit to 
the LASS data [26] . The two sets of values are consistent 
within errors and lead to similar results. 

To account for the three-body kinematics in the anal- 
ysis of B'^ — > (pKir decays, we multiply the amplitude 
squared |A/j(m)p by the phase-space factor F{m): 



F{m) = 2 X m X [toLx(™) " '"min(™)] . (26) 



where 



^max cind TOujjjj 



are the maximum and minimum 



values of the Dalitz plot range of {m^K is the ipK 

invariant mass) at any given value of mxTn see kinematics 
section of Ref. [l^ . Due to slow dependence of the factor 
in Eq. (|26p on m in any small range of m, the difference 
of this approach from the quasi-two-body approximation 
is small. 



C. Interference effects 

The differential decay width discussed in Eq. ([7]) in- 
volves interference terms between resonances with differ- 
ent spins J. These interference terms have unique angu- 
lar and mass dependences which cannot be factorized in 
the full distribution. We can parameterize the mass and 
angular amplitude for each spin state J as follows: 
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TABLE III: Parameterization of the angular distribution in Eq. (|8} in the B" — > (p{Ktt)j decays where three resonance final 
states with spin J = 0, 1,2 are considered. The common constant is quoted for each decay mode and is omitted from each 
individual function below. The three helicity angle parameters {Ti.i, Tt2, ^) are discussed in the text. 







(pK* (892)" 


V3A'2*(1430)° 




J — U 


r — 1 


7 — 9 


common constant 


3/47r 


9/S-K 


15/327r 


fij (Hi.Ha, $) 




rti ri2 


- l)^Hi 


/2j (Hi.Ha, $) 





i(l-H?)(l-Hi) 


3H?(l-H?)(l-'Hi) 


/3J {-Hi, H2, $) 





i(l-H?)(l-Hi)cos2$ 


3H?(1 -Hi)cos2$ 


/4j {-Hi, H2, <&) 





-i(l-H?)(l--Wi)sin2$ 


-6H?(1 -'Wi)sin2$ 


/5,/(Hi, n2, $) 





ViUi^i - nln2^/i - nl cos$ 


Vl - (3-W^ - 1) n2 x/1 - cos $ 


/6,/(Hl, H2, -l-) 





-V2nly/l-njn2^/l-nl sin* 


- V6Hi - 'Hi (Sni - 1) - sin $ 




FIG. 4: Intensity |Mj(mif7r)|'^ (a) and phase arg(Mj(mKir)) (b) of the invariant amplitudes for J = (solid), J — 1 (dashed), 
and J = 2 (long-dashed) Kn contributions as a function of the invariant Kn mass ttiktv The taller two arrows indicate the 
low ttia'tt region, while the shorter two arrows indicate the high mjfir region. The relative intensity of the amplitudes is taken 
from Fig. [2] while the absolute intensity is shown in arbitrary units. 



Ao(m;f,,0i,02,$) = Y°{Hi,<i>)Y^"{-n2,0)Mo{mK.)Aoo (27) 
Ai{mK.,0i,e2,<i>) = J2 Y^^{ni,^)Y^-\-n2,0)Mi{mK^)Aix (28) 

A=0,±1 

A2{mK.,0i,e2,<i>) = Y^\ni,'i>)Y^-\~n2,0)M2{mKn)A2x. (29) 

A=0,±1 



The interference will appear in the angular-mass distri- two interference terms, properly normalized, are defined 
buttons as 2'Sie{Ai{mK7T , 0i,92,^)A*{mKTr, 6*1, 6*2, $)). As for J = 1 and J = 2: 
we can see from Fig. 21 the overlap between the P- and 
D-wave Ktt contributions is negligibly small, and we will 

consider only the interference between the J ~ and g 

J = 1, or J = and J 2 amplitudes. The resulting jMe{AjA^ ^ y. ^-^^^^^^ /,j(Hi, H2, $) , (30) 



vmjii^iAo 
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TABLE IV: Parameterization of the Ktt invariant mass distri- 
bution in the B'' — > Lp[K-K)j decays where three resonance fi- 
nal states with spin J = 0, 1, 2 are considered. The resonance 
mass mj, width Fj [l^, [23, [2^, interaction radius r, scat- 
tering length a, and effective range b are considered [25l. [2^. 
Combined errors are quoted, except for {K-k)q'^ where the sys- 
tematic errors are quoted last while the central values and sta- 
tistical errors have been updated [26[ with respect to Ref. . 







if* (892)° 


7^2* (1430)° 




J = 


J= 1 


J = 2 


mj (MeV) 


1435 ± 5 ± 5 


896.00 ± 0.25 


1432.4 ± 1.3 


Tj (MeV) 


279 ± 6 ± 21 


50.3 ±0.6 


109 ±5 


r (GeV-i) 




3.4 ±0.7 


2.7 ± 1.3 



a (GeV"^) 1.95 ±0.09 ±0.06 
6 (GeV~^) 1.76 ±0.36 ±0.67 



afj(mK^) = 
agjirriKT,) = 



The above terms are specific to the i?" decays and are 
denoted with superscript "— " , as introduced in Eqs. ((T]- 
[S]). The interference distributions for the B° decays are 
described by the same Eq. pO|) , but replacing A by A and 
a~j by a^j and with definitions given by Eqs. ([3TH33|l . 
replacing superscript "— " by "+" . 

The main difference now is that the a~j{mKTT) parame- 
ters, as defined for i — 7, 8, 9, have a different dependence 
on mass to those defined for i — 1-6 in Eqs. (Pl [T?|l . This 
dependence now includes the phase of the resonance am- 
phtude as a function of mass. This dependence becomes 
crucial in resolving the phase ambiguities. 

As can be seen from Eq. ([5]) and Eqs. ([MH)). for 
any given set of values (^y j, A0|| j, A0^j) a sim- 
ple transformation of phases, for example (27r — j, tt — 
— A(/)|| J, — A(/)_L,/), gives rise to another set of val- 
ues that satisfy the above equations in an identical man- 
ner. This results in a four- fold ambiguity (two- fold for 
each of and decays). At any given value of 
ttikt: the distributions, including the interference terms 
in Eqs. ([5^ and ([M)) . are still invariant under the 
above transformations if we flip the sign of the phase, 
arg(Mj(mif^)AfQ (TOx7r)e~"'o-'). At a given value of 
this phase has to be determined from the data and we 
cannot resolve the ambiguity. However, the mass depen- 
dence of this phase is unique, given that the parameters 



where the angular dependence is defined in Table fVl and 
a~j{mK-w) are defined for z = 7, 8, 9 as: 



(31) 
(32) 
(33) 

I 

5qj are constant. Therefore, the two ambiguous solutions 
for each B^ and decay can be fully resolved from the 
niK-jT dependence of the angular distributions in Eq. (j30p . 

This technique of resolving the two ambiguous solu- 
tions in _B — > VV decays has been introduced in the 
analysis of — + J/ipK*'^ decays 27] and has been used 
in BaBAR's earlier analysis of both B ipK*° and (pK*"^ 
decav s [20l . l2lj . It is based on Wigner's causality prin- 
ciple (28|, where the phase of a resonant amplitude in- 
creases with increasing invariant mass, see Eq. p3|) . As 
a result, both the P-wave and Z?-wave resonance phase 
shifts increase rapidly in the vicinity of the resonance, 
while the corresponding 5-wave increases only gradually, 
as seen in Fig. S) 



D. Time-dependent distributions 

Measurement of the time-dependent CP asymmetry 
A{At) in the decay of a neutral B meson to a CP eigen- 
state, dominated by the tree-level b c amplitude or 
by the penguin b ^ s amplitude, such as i3° {cc)Kg 
or B° — > {ss)Kg, where (cc) and (ss) are charmonium or 
quarkonium states respectively, gives an approximation, 
/3off, to the CKM Unitarity Triangle angle (3 [2^]. The 



7^ !iie{Mj{mK.)M*imK^)e-'^oj) 
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TABLE V: Parameterization of the angular distribution in Eq. (|30p . Interference between either J = and J = 1, or J = 
and J = 2, contributions in the tp{K-K)j decays is considered. The common constant is quoted for each decay mode and 

is omitted from each individual function below. The three helicity angle parameters (7ii, 7^2, "!>) are discussed in the text. 











J = 1 


J = 2 


common constant 








HiUl 


{3hI - i)nl 


fsj (Hi, Ha,*) 


^^l-Ul^l-HiU2 cos$ 


V6a/1 - H? T-ti \/l - T-tl ^2 cos $ 


,f9j(Hl, Hi, <&) 


-^^l-nl^/l-HiHi sin-l- 


- V6\/l - T-i-l Hi - nl Hi sin $ 



CP asymmetry is defined by 



A{At) 



N{At,B^,^)-NiAt,B^,,^) 
N{At,Bl^) + NiAt,B°,,^) 

= S'sin(AmBAi) - Ccos(AmBAi) , (34) 



and 



- sin(2/3cff) 



— \sm — 7 



TJCP X 5/VT 





9^ 


)/ 





C2, 



(35) 



where N{At,B^^^) or iV(At,B^ag) is the number of 
events observed to decay at time At, in which the fla- 
vor of the B meson opposite to that of the decaying B at 
At = (referred to as the flavor "tag") is known to be 
or B^ respectively, r/cp = ±1 is the CP eigenvalue of 
the final state; amplitudes A and A describe the direct 
decays of 5° and B'^ respectively to the final state; and 
AniB is the mixing frequency due to the difference in 
masses between the B meson eigenstates. We use a con- 
vention with A — rjcp X A in the absence of CP violation. 
The above asymmetry follows from the time evolution of 
each flavor: 



N{At,Bl 



-|At|/TB 



X (1 + 5sin(Ami3At) - C cos(Ams Ai)) (36) 



iV(At,i?° )cx 



j-|At|/rB 

Atb 



X (1 - S'sin(AmBAt) + Ccos(ATOBAt)) . (37) 

The B^-B^ mixing parameters q and p can be ex- 
pressed to a good approximation using the Wolfenstein 
phase convention within the Standard Model [3| '■ 



(38) 
(39) 



arg ( ^ ) = -2/3 



= 1. 



The value of sin 2/3 from the charmonium b ^ c decays, 
which defines the phase of the mixing diagram in the 
Wolfenstein parameterization, is well measured [sl. [Tst: 



sin 2/3 ==0.681 ±0.025 
or 2/3 = (0.75 ± 0.03) rad 



(40) 



where the phase ambiguity of /3 in the range [0, tt] has 
been resolved using vector-vector charmonium B de- 
cays [13 and the decay B" ^ K+R-R^ [s^. Should 
there be a new physics contribution to the mixing dia- 
gram, its effect is absorbed into the definition of [3 in 
Eq. ([55]) , which should be valid for the analysis discussed 
in this paper. New physics effects in the 5 — > c ampli- 
tude are unlikely to be significant as this transition is not 
suppressed in the Standard Model. Therefore the com- 
parison of sin 2/3 in Eq. ([3D]) with sin 2/3eff measured in 
b s transitions would be a test of new physics in the 
penguin B decays. 

There is an alternative sign convention for the choice 
of the direct-CP violation parameter defined by: 



C 



Acp = 



|A|2 + |A|2- 



(41) 



Given the approximation in Eq. p9p and our phase 
convention in Eq. (|38p . the value of 5 can be expressed 
as: 



5 = ^/1 



-2/3 + arg ( - 



(42) 



Therefore, we have S = —rjcp x sin 2/3 when A = 
rjcp X A. When we measure S in Eq. (|42)) with b ^ s de- 
cays, we can safely assume that the value of /3 has been 
measured in charmonium decays, as given in Eq. (|40p . 
Therefore, we are ultimately interested in the measure- 
ment of arg{A/A). Any large deviation from arg(?7cp) 
would be a signal of new physics. 

In the study of the time-evolution in Eq. ([M]) we can 
use the decay B^ — > ip{RgT:^)Q^ with an S-wave Rtt 
contribution. This final state is a CP-eigenstate with 
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Vcp = +1 as we discuss below. However, the situation 
is more complicated in the general case of (pKgTr'-' 
decays where the final state is no longer a CP-eigenstate. 
The amplitude for this decay is a superposition of CP 
eigenstates. 

In Ref. 31] it was shown that the CP quantum num- 
bers are independent of the K^tt'^ system, that is inde- 
pendent of J, for the decay i?'^ {cc)KgTr'^. The same 
analysis applies to {ss)KgTr'^. The CP parity is 

defined only by the (ss) spin alignment A (an alternative 
analysis that introduces the eigenstate of the transversity 
T is sometimes used to separate the CP eigenstates [Slj). 
For example, all longitudinal decays — > (pKgTr^, cor- 
responding to A = 0, are CP-even, including the decay 
93(^^5 7r'^)3°. Overall, we conclude that in the decay 
B" ip{KgTT^)j we have three amplitudes with definite 
CP: 









Vcp = - 


fl 


(43) 
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= {Aj+ 


+ Aj_)/V2 


Vcp = - 


fl 


(44) 




= iAj+ 




Vcp = - 


-1. 


(45) 



Similarly, ijcp = — 1 for the B decay to the final state 
foK gTr'^, which will be considered as a background decay 
in our analysis. We do not discuss the CP properties 
of the interference terms with the product of amplitudes 
of different CP; these terms are integrated over in our 
analysis of the time evolution. These terms could be 
considered in a future experiment with higher statistics. 

We can express the time-evolution coefficient in the 
decay B° ip{K^7r°)^°, as 



Soo = -^Jl-Alo^ sin(2/3 + 2A0oo) , (46) 



where ^00 — Acpo- For the decays B'^ (p{Kg7r'^) with 
if* (892)0 (J ^ K*{U30)° (J = 2) intermediate 

resonances there are three time-evolution terms, one for 
each amplitude, if we ignore the interference terms: 



Sjo 
Sj± 



1 - Ajo X sin(2/3 -I- 2A(5oj + 2A0O 



1 - A^j,. X sin(2/3 -I- 2A(5oj - 2A(/)|| j -t- 2A0oo) 



1 - A^j^ X sin(2/3 + 2A(5o.7 - 2A<j)±j + 2A0oo). 



(47) 
(48) 
(49) 



The three corresponding direct-CP violation terms Ajo, Aji\, and Aj±_ can be obtained from the direct-CP violation 



and polarization parameters measured in the i? ^ — > (p(K^'K^) decays: 

Acpj + Aqpj 



A.JO = 
Aj± = 

M = 



1 + Acpj x A°cpj 

Acpj + A^pj 
1 + Acpj x A^pj 

Acpj - Jlj x (Acpj + A°cpj) - f±j x {Acpj + A^pj) 
1 - /lj X (1 + Acpj x A°cpj) - f±j x (1 + Acpj x A^pj) ' 



(50) 
(51) 
(52) 



r 



As can be seen from Eqs. (|46ff52| . there is only one 

parameter A(/)oo that is not measurable in the B ^ 
(p{K^TT^) decays. Therefore, the above parameterization 
allows us to measure A0oo from the time evolution of the 
(p{KgTT'^) decays while all other parameters are 

measured in the mode _B ^ — > ip{K^TT^), which has a 
significantly larger reconstructed yield. 

The angular distributions in Eq. ([7]) can be simplified 
after integrating over the angle The resulting angular 
distribution will not have interference terms between dif- 
ferent amplitudes for a given J. This makes the time evo- 
lution parameterization relatively simple with just two 



terms: longitudinal {/lj) and transverse (1 — /lj) polar- 
ization. The longitudinal time-evolution is parameterized 
by the Sjo coefficient, and the transverse time-evolution 
is parameterized by the expression 



'JT 



1 - 



^ f±j + Acpj x A^pj) 
~ 1 - /lj X (1 + Acpj x A^^pj) 

f±j X (1 + Acpj x A^pj) 
1 - /lj X (1 + Acpj x A°cpj) 



X S. 



jj- 



X5;||. (53) 



In a similar manner, the longitudinal direct CP-violation 
term is parameterized by the Cjo = —Ajo coefficient, 
and the transverse direct CP-violation term is parame- 



15 



terized by 



^ f±J X (1 + -^cpj X A^pj) 



1 - Jlj X (1 + AcPJ X A%j) 
f^j X (1 + ylppj X A^pj) 
1 - /l.7 X (1 + AcPJ X y^O,pj) 



X ^ 



(54) 



As an example, let us consider the J = 1 case. Since 
Ai± and have opposite CP-parity and it has been 
measured that /^i ~ /yi = (1-/±i-/li) BEIIS, in 
the Standard Model we expect to a good approximation 
Sit ^ 0, S'lo = sin 2/3, and Cit = Cm = 0. 



III. EVENT RECONSTRUCTION 

We use a sample of (465.0 ± 5.1) million e+e^ 
r{4S)-^ BB events collected with the BaBAR detec- 
tor [34I at the PEP-II e+e" asymmetric-energy storage 
rings. The center-of-mass system of the T{AS) resonance 
is boosted providing roughly 250 /im average separation 
between the two B meson decay vertices. The e^e~ 
center-of-mass energy ^/s is equal to 10.58 GeV, corre- 
sponding to the T(45') resonance. 

We fully reconstruct the (^(1020) 

[K^ K^){ K -k) candidates with two {Kn) final states, 
KliT° and K^t:' 

are reconstructed in the final states 



The neutral pseudoscalar mesons 
TT+TT" and 

TT*^ — !■ 77. The dominant background in our analysis 
comes from e'^e'^ — > qq production {q = u,d,s,c). A 
data sample equivalent in luminosity to 12% of the on- 
T(4S')-resonance sample has been collected with y/s 40 
MeV below the T(4S') resonance (off-resonance data) for 
studies of this background. A detailed GEANT4-based 
Monte Carlo (MC) simulation [33] of the detector has 
been used to model all processes. This simulation has 
been extensively tested and tuned with high-statistics 
validation samples. 

Momenta of charged particles are measured in a track- 
ing system consisting of a silicon vertex tracker (SVT) 
with five double-sided layers and a 40-layer drift cham- 
ber (DCH), both within the 1.5-T magnetic field of a 
solenoid. Identification of charged particles (PID) is pro- 
vided by measurements of the energy loss in the tracking 
devices (dE/dx) and by a ring-imaging Cherenkov detec- 
tor (DIRC). Photons are detected by a CsI(Tl) electro- 
magnetic calorimeter (EMC). We use minimal informa- 
tion from the muon identification system (IFR) to make 
a loose veto of the charged muon tracks. 

We require all charged-particle tracks (except for those 
from the Kg —> tt+tt" decay) used in reconstructing the 
B candidate to originate from within 1.5 cm in the x-y 
plane and 10 cm in the z direction from the nominal beam 
spot. We veto leptons from our charged-particle track 
samples by demanding that tracks have DIRC, EMC, and 
IFR signatures that are inconsistent with either electrons 



or muons. Further pion and kaon PID requirements are 
based on a likelihood selection developed from dE/dx 
and Cherenkov angle information from the tracking de- 
tectors and DIRC, respectively. The typical efficiency of 
PID requirements is greater than 95% for charged tracks 
in our final states. Photons are reconstructed from en- 
ergy deposits in the electromagnetic calorimeter that are 
not associated with a charged track. We require that all 
photon candidates have an energy greater than 30 MeV 
in the EMC. 

The invariant mass of the candidate Kg is required 
to lie within the range |m^+^- — mKo\ < 12 MeV. We 
perform a vertex-constrained fit to require that the two 
tracks originate from a common vertex, and require that 
the lifetime significance of the Kg be r/ar > 5, where 
T and a-r are the Kg lifetime and its uncertainty deter- 
mined from the vertex-constrained fit. For the Kg candi- 
dates, we also require the cosine of the angle between the 
flight direction from the interaction point and momen- 
tum direction to be greater than 0.995. The efficiency of 
Kg selection requirements is about 90%. 

We select neutral-pion candidates from two photon 
clusters with the requirement that the 77 invariant mass 
satisfy 120 < m^j < 150 MeV. The mass of a tt" can- 
didate meetin g th is criterion is then constrained to the 
nominal value [18] and, when combined with other tracks 
or neutrals to form a B candidate, to originate from the 
B candidate vertex. This procedure improves the mass 
and energy resolution of the parent particle. The purity 
of A'O and tt" selection is 92% (88%) and 90% (68%), 
respectively, in the signal sample (combinatorial back- 
ground) based on MC simulation studies. 

We identify B meson candidates using two main kine- 
matic variables, beam energy-substituted mass m-Es- 



TOES = [(s/2 + pr • PbY/Et - pi] 



2.11/2 



(55) 



and the energy difference AE: 

AE = {EtEb VrVB s/2)/^s, (56) 

where (i?B,Ps) is the four- momentum of the B can- 
didate, and (i^rjPr) is the e+e^ initial-state four- 
momentum, both in the laboratory frame. Both vari- 
ables are illustrated in Fig. O The distribution of AE is 
expected to peak at zero and toes at the B mass around 
5.280 GeV. The AE resolution is dominated by the de- 
cay product energy and momentum measurements and 
is typically 34 and 20 MeV for the subchannels with and 
without a tt" , respectively. The typical toes resolution is 
2.6 MeV and is dominated by the beam energy uncertain- 
ties. We require toes > 5.25 GeV and \AE\ < 100 MeV 
to retain sidebands for later fitting of parameters to de- 
scribe the backgrounds. 

The requirements on the invariant masses of the res- 
onances are 0.99 < to^^ < 1.05 GeV and 0.75 < 
TOj^o^o < 1.55 GeV, 0.75 < to^±^t < 1.05 GeV, or 

1.13 < mx±TrT < 1-53 GeV for the ip and K*^, respec- 
tively. Here we separate the K^tt^ invariant mass into 



two ranges for later fitting. The two ranges simplify the 
fit configuration and allow us to test the nonresonant 
B ipKn contribution independently, therefore provid- 
ing an independent crosscheck of the parameterization. 

To reject the dominant e'^e" —^ qq background, we 
use variables calculated in the center-of-mass frame. We 
require |cos0t| < 0.8, where Ot is the angle between 
the i3-candidate thrust axis and that of the rest of the 
event. The angle 9t is the most powerful of the event 
shape variables we employ. The distribution of | cos 6t \ 
is sharply peaked near 1 for combinations drawn from 
jet-like qq pairs and is nearly uniform for the isotropic 
B-meson decays. Further use of the event topology is 
made via the construction of a Fisher discriminant 
which is subsequently used as a discriminating variable 
in the likelihood fit. 

Our Fisher discriminant is an optimized linear combi- 
nation of the remaining event shape information, exclud- 
ing cos^T- The variables entering the Fisher discrimi- 
nant are the angles with respect to the beam axis of the 
B momentum and B thrust axis, and the zeroth and sec- 
ond angular moments Lo,2 of the energy flow about the 
B thrust axis, all calculated in the T{AS) center-of-mass 
frame. The moments are defined by 

L, =Y,P^x\cos9,W (57) 

i 

where di is the angle with respect to the B thrust axis 
of track or neutral cluster i, Pi is its momentum, and the 
sum excludes the B candidate. The coefficients used to 
combine these variables are chosen to maximize the sep- 
aration (difference of means divided by quadrature sum 
of errors) between the signal and continuum background 
distributions of Lj, and are determined from studies of 
signal MC and off-peak data. We have studied the opti- 
mization of J- for a variety of signal modes, and find that 
the optimal sets of coefficients are nearly identical for 
all. Because the information contained in is correlated 
with I cos 9t I , the separation between signal and back- 
ground is dependent on the | cos 9t \ requirement made 
prior to the formation of !F. The !F variable is illustrated 
in Fig. O 

In order to establish that the MC simulation repro- 
duces the kinematic observables in the data, such as 
Toes, ^E, and JF, we use high-statistics i3°-meson de- 
cays with similar kinematics and topology. For exam- 
ple, in Fig. [5]we illustrate reconstructed 5° D~7r+ 
(if+7r^7r^)(7r+) decays. There is good agreement be- 
tween data and MC. The deviations in the means of the 
distributions are about 0.7 MeV for toes, 5 MeV for AE, 
and negligible for We take these corrections into ac- 
count when we study the B° ipKn decays. 

The B background contribution is generally found to 
be small due to selection on the narrow tp resonance, 
PID requirements on the kaons, and good momentum 
resolution, important in particular for AE. We remove 

TT^ signal candidates that have decay prod- 
ucts with invariant masses within 12 MeV of the nominal 
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FIG. 5: Validation of kinematic variables with the high statis- 
tics sample of B° — > D~tt'^ — > {K'^TT~Tv~){Tr'^) decays. Pro- 
jections onto the variables mES, and J- are shown from 
top to bottom for MC (a), (b), (c), and data (d), (e), (f). 
The dots with error bars represent the MC simulation (left) 
or data (right). The long dashed lines represent the signal 
and the solid lines show signal-plus-background parameteriza- 
tion. A small fraction of combinatorial background is present 
(dashed lines) due to combinations of other B decays in the 
MC and also qq continuum in the data. 



mass values for Df or D^. This removes the background 
from the B DfK, D^K, D±7r, and D^tt decays. To 
reduce combinatorial background in the (pKgir'^ 
analysis with low-momentum tt^ candidates, we require 
Hi < 0.8. Certain types of B background, such as po- 
tential B° fo{980)KTT, cannot be distinguished on an 
event-by-event basis, and so cannot be removed by ve- 
toes. We incorporate these contributions into the fit. 
The remaining B background events were found to be 
random combinations of tracks and can be treated as 
combinatorial background, similar to random tracks from 
qq production. 

When more than one candidate is reconstructed, which 
happens in 5% of (pK^n^ and 10% of (fKgTT° events, we 
select one candidate per event based on the lowest value 

of the of the four-track vertex for B ^ ^ ipK^n^ or 
of the fitted s''' ipK'^Tr^ decay tree. 

In the self-tagging i?-decay mode B^ ^ ipK^ii^ , we 
define the 6-quark fiavor sign Q to be opposite to the 
charge of the kaon candidate. For each reconstructed 

B ^ ^ LpK^Ti'^ signal candidate (-Bgig) we use the re- 
maining tracks in the event to determine the decay vertex 
position and flavor of the other B decay, i?tag- A neural 
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network based on kinematic and particle identification in- 
formation assigns each event to one of seven mutually ex- 
clusive tagging categories (ctag) [s^l, including a category 
for events in which a tag flavor is not determined. The 
_B-flavor-tagging algorithm is trained to identify primary 
leptons, kaons, soft pions, and high-momentum charged 
particles from the other B and correlate this informa- 
tion to the B flavor. The performance of this algorithm 
is evaluated using a data sample (Bflav sample) of fully- 
reconstructed B^ D^*^~iT~^ /p'^ /ai decays. The effec- 
tive tagging efficiency is measured to be (31.2 ± 0.3)%. 

We determine the proper time difference, At, between 
Bsig and Btag from the spatial separation between their 
decay vertices. The Btag vertex is reconstructed from the 
remaining charged tracks in the event, and its uncertainty 
dominates the At resolution a At- The average proper 
time resolution is (cta*) ~ 0.7 ps. Only events that satisfy 
|At| < 15 ps and 0.1 < cta* < 2.5 ps are retained. 

Overall, selection requirements discussed here have 
been optimized to retain large signals and wide side- 
bands of observables for later fitting, as discussed in 
the next section. Statistical precision of the measure- 
ments was the main optimization factor, while individ- 
ual systematic uncertainties were kept small compared 
to statistical errors. After applying all selection crite- 
ria and using the measured branching fractions summa- 
rized in Table HI we expect to observe about 177 (15) 

W° ip{ '^V)*o, 473 (34) W° if '^'*(892)", and 156 

(9) 's'" ^ ip ^'^(1430)" events in the (pK^n^ (v'^s'^") 
channel. The larger reconstruction efficiency and sec- 
ondary branching fractions in the ipK^n^ channel result 
in the dominance of this decay mode in the signal param- 
eter measurements, except for the measurement of the 
A0OO parameter, which is possible only with the (fiKgir^ 
channel. 



IV. ANALYSIS METHOD 

We use an unbinned, extended maximum-likelihood 
(ML) fit [ij to extract the 27 parameters defined in Ta- 
ble |TT1 which describe three decay channels (12 in either 
B° (^i^* (892)0 ipK;{U30)°, and three in ip{KTT)*° 
decays). We perform a joint fit to the data for three 
modes: (pKgTT^; (fK^n^ in the lower Ktt mass range 
(0.75 - 1.05 GeV); and (pK^n^ in the higher Ktt mass 
range (1.13 - 1.53 GeV). To simplify treatment of the 
likelihood function, we separate the two ranges of the 
K^TT^ invariant mass. However, in the joint fit the like- 
lihood function C is written as a product of three inde- 
pendent likelihood functions, one for each of the above 
three modes, as discussed below. 

Due to the relatively low statistics, we simplify the 
angular analysis in the -8° ipKgiT^ decay mode. We 
integrate the angular distributions in Eq. ([7]) over the 
angle $. Therefore only the longitudinal polarization 
fractions /lj, the yields, and the time-evolution are mea- 



sured with the i?" (pKgTT^ decays. We constrain the 
relative signal yields for the same spin-J contributions 
in the two subchannels ip{K^n^)j and ip{KgTr'^)j taking 
into account the isospin relationship, daughter branching 
fractions, and reconstruction efficiency corrections. The 
isospin relationship requires that the K**^ K^n^ and 
K^tt'^ fractions are 2/3 and 1/3 of the total Kn decay 
rate, and we ignore any isospin violation as being neg- 
ligible for the measurements in this analysis. All other 
signal parameters in Table [IT] are constrained to be the 
same when they appear in both channels. 

A. Likelihood function 

The likelihood function for (pKgn'^ is written as: 

£ = n exp i-N,) n [ J2 n,f^V]{x,; C, ^) ) , (58) 

where rij is the unconstrained (except if noted otherwise) 
number of events for each event type j, /J is the fraction 
of events of component j for each tagging category c, 
Nc ~ "-J the number of events found by the fit 

for tagging category c, and Vj{xi] the probability 
density function (PDF). 

The data model has five event types j: the signal 
5° (f{KTr)j with J = 0,l,2, a possible background 
from i?" fo{980)K*'^, and combinatorial background. 
The combinatorial background PDF is found to ac- 
count well for both the dominant light qq events and 
the random tracks from B decays. Each event candi- 
date i is characterized by a set of 10 observables Xi — 

{toes, AE, JT, rriKn, rUj^x, "Hi, "^2, Ctag, At, CTAtji, 

the kinematic observables rriES, AE, T , the K* and <p 
invariant masses raKix and rrij^j^, the helicity angles Hi, 
Ti.2, the flavor tag Ctag, the proper time difference At and 
its event-by-event error a At ■ The PDFs are split into the 
seven tagging categories. The polarization parameters 
quoted in Table are denoted by and the remain- 
ing parameters by ^. Most of the ^ parameters, except 
for A0OO, appear also in the likelihood function for the 

i? — > tpK^TT^ decays. 

The likelihood function for _B ^ — > ipK^n^ decays, in 
either the lower or the higher Kn mass range, is 

(59) 

where the index j represents the three event types used in 
our data model: the signal B'^ (p{Ktt)'^ (j = 1) which 
combines the two dominant modes in a given mass range 
(V?^* (892)0 ip{KTT)f in the lower or ipK^ {1A30)° 
and (p{Ktt)q^ in the higher mass range), a possible back- 
ground from B° — > /o(980)X*o f^j _ 2), and combinato- 
rial background (j = 3). The superscript k corresponds 
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to the value of the flavor sign Q = ±1 and allows for 
a CP- violating difference between the B° and decay 
amplitudes {A and A). 

In the signal event type, the yield and asymmetry, risig 

and Acp, of the 'b^° (p{K^Tr^)j mode with J = 1 
in the lower mass range or J = 2 in the higher mass 

range, and those of the B ^ ^ lp{K^tt'^)'^ mode are 
parameterized by applying the fraction fj}' of the (p{Ktt)j 
yield to nj. Hence, rigig = x /i"*" -I- nj~ x Acp ~ 
(rij^ X /i+ — X ^^)/nsig, and the ip{K'^n^)Q^ yield 
is n'l X {1 — ) + Ui X {1 — This treatment is 

necessary to include interference between the two decay 
modes as we discuss below, while we ignore interference in 
the i?" — > ifiKgir'^ channel due to low statistics. The PDF 
is formed from the following set of observables x.^ ={toeSj 
AE, T , rriK-w, t^^kk^ ^2, Q}, and the dependence 
on /i*^ and polarization parameters ( = {/lj, f±j, 
'/>±j, Sqj, Al;pj, Acpj, A(?!)||j, A(j)±j, A6oj} is relevant 
only for the signal PDF Vi- 

The remaining PDF parameters ^, in both the B ^ 
(pK^-K^ and ipK^ir^ channels, are left free to vary in 
the fit for the combinatorial background and are fixed to 
the values extracted from MC simulation and calibration 
B — > -Dtt decays for the other event types. We minimize 
the -21n£ function using MINUIT ^ in the ROOT 
framework [36| . The statistical error on a parameter is 
given by its change when the quantity — 21n£ increases 
by one unit. The statistical significance is taken as the 
square root of the difference between the value of —2 In £ 
for zero signal and the value at its minimum. We have 
tested this procedure with simulated samples and found 
good agreement with the statistical expectations. 



B. PDF parameterization 

The PDF "P^Xx,; C; for a given i?'" •pK^-K^ 
or (pK'giT^ candidate i is taken to be a joint PDF for 
the helicity angles, resonance masses, and Q, and the 
product of the PDFs for each of the remaining variables. 
The assumption of negligible correlations in the selected 
data sample among the discriminating variables, except 
for resonance masses and helicity angles where relevant, 
has been validated by evaluating the correlation coeffi- 
cients. This assumption was further tested with the MC 
simulation. 

For the parameterization of the signal PDFs we use 
double-Gaussian functions for signal AE and mES- For 
the background we use low-degree polynomials as re- 
quired by the data or, in the case of mESi an empirical 
threshold ARGUS function [13 : 

f{x) cx xy/l-x"^ exp[-Ci(l - x^)] , (60) 

where x = WEs/^'beam and is a parameter that is 
determined from the fit with a typical value of about 25. 



For both signal and background, the Fisher distribu- 
tion T is described well by a Gaussian function with dif- 
ferent widths to the left and right of the mean. For the 
combinatorial background distribution, we also include a 
second Gaussian function with a larger width to account 
for a small tail in the signal J- region. This additional 
component of the PDF is important, because it prevents 
the background probability from becoming too small in 
the region where signal lies. 

A relativistic spin- J B-W amplitude parameterization 
is used for the resonance masses 0, [M] , except for the 
{Ktt)^ rriKTr amplitude, which is parameterized with 
the LASS function [H, Hi]. The latter includes the 
i^Q(1430)° resonance together with a nonresonant com- 
ponent. The detailed treatment of the invariant mass 
distribution is discussed in Section HlBI We found that 
no additional correction to the Ktt invariant mass param- 
eterization is necessary because resolution effects of only 
a few MeV are negligibly small compared with the reso- 
nance widths. On the other hand, we convolve resolution 
effects into the rrij^j^ parameterization. 

The background parameterizations for candidate 
masses include resonant components to account for res- 
onance production in the background. The background 
shape for the helicity parameterization is also separated 
into contributions from combinatorial background and 
from real mesons, both fit by low-degree polynomials. 

The mass-helicity PDF is the ideal distribution from 
Eqs. (|27H30p . multiplied by an empirically-determined 
acceptance function Q{Hi,H2,^) = GiiTi-i) x 02(7^2), 
which is a parameterization of the relative reconstruc- 
tion efficiency as a function of helicity angles. It was 
found with detailed MC simulation that resolution ef- 
fects in the helicity angles introduce negligible effects in 
the PDF parameterization and fit performance, and they 
are therefore ignored. The angles between the final state 
particles and their parent resonances are related to their 
momenta. The signal acceptance effects parameterized 
with the function G{Ti.i, Ti.2, $) are due to kinematic cor- 
relations, whereas the detector geometry correlations are 
negligible. Therefore the above uncorrelated parameter- 
ization as a function of two helicity angles was found to 
be appropriate and was validated with detailed MC sim- 
ulation. 

Momentum in the laboratory frame is strongly cor- 
related with detection efficiency. Thus we have accep- 
tance effects in the helicity observables Hi , most evident 
for the large values of Tii corresponding to the slow tt 
from the K* meson. However, these acceptance effects 
are not present for the $ angle; there is no correlation 
with the actual direction with respect to the detector, 
which is random for the B decays. The acceptance ef- 
fects for the two helicity angles Hi and H2 are shown in 
Fig. [SI We obtain the acceptance functions from the fit 
to the signal MC helicity distribution, with the known 
relative components of longitudinal and transverse am- 
plitudes generated with B° (fK*{892)'^ MC. The D^^- 
meson veto causes the sharp acceptance dips around 0.8 
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FIG. 6: Angular acceptance functions for Tii in ipK^-rv^ (a), in (^A'gTr" (b) , and for 7^2 (c). These plots show only relative 
efficiency between different helicity points with arbitrary y-axis units. The D^j-meson veto causes the sharp acceptance dips 
near Tii — 0.8 seen in (a). 



in the Gii'Hi) function in the B ^ ^ ipK^n^ analysis. 

The interference between the J = 1 or J = 2 and 
the S'-wave (Ktt) contributions is modeled with the term 
2^e{AjAQ) in Eq. ([30]) with the four-dimensional an- 
gular and rriK-K dependence, as discussed in detail in 
Section III CI It has been shown in the decays i?° 
J/^{Ktt)1^ and B± tt^{Kit)1^ ^ that the ampli- 
tude behavior as a function of m^vr is consistent with 
that observed by LASS except for a constant phase shift. 
Integrating the probability distribution over (Tii, 7^2, <&), 
the interference term 2^e{AjAQ) should vanish. How- 
ever, as we introduce the detector acceptance effects 
on (7^1,7^2), the interference contribution becomes non- 
zero. The total yield of the two modes is corrected before 
calculating the branching fractions. The effect can be es- 
timated by comparing the integral of (/?i4r|(1430)'', 
_B" — > lp{Ktt)'^ and the interference probability contri- 
bution. We find that the interference term accounts for 
3.5% of the total yield. Accordingly, we scale the yields of 
B° 93/4:2* (1430)° and B" ^ v{Kt:)1^ modes by 96.5% 
while calculating the branching fraction. This effect is 
neghgible for the B° ipK*{%%2f decay due to the rel- 
atively small fraction of the B^ lp{K'k)^ contribution 
to the lower Kt: mass range. 

The parameterization of the nonresonant signal-like 
contribution -> f^K*^ {K+K-)K*" is identical 



to the signal in the primary kinematic observables tties, 
A_E, J^, and the rriK-K mass but is different in the angu- 
lar and rrij^j^ distributions. For B'^ /qA'*", the ideal 
angular distribution is uniform in $ and 7^2 and is pro- 
portional to Ti-i , due to angular momentum conservation. 
We use a coupled-channel B-W function to model the 
K^K~ mass distribution for the /o f38l|. The broad in- 
variant mass distribution of /□ compared to the narrow Lp 
resonance was found to account well for any broad ttIj^j^ 
contribution. This PDF parameterization is further var- 
ied as part of the systematic uncertainty studies. The 
mjiTT distribution in JqK*^ is parameterized as a 

J = 1 contribution in the lower mass range and as a J = 
contribution in the higher mass range, see Section [II Bl 

For the B'^ — > ipK^ii'^ mode, the $ angle is integrated 
over, so that no interference terms appear in the fit. 
An additional PDF for the Ai distribution is used for 
both the signal and background, which is discussed next. 
The treatment of other observables is similar to those of 

Time-dependent CP asymmetries are determined using 
the difference of B^ meson proper decay times At = ^sig — 
itag, where igig is the proper decay time of the signal B 
(Bsig) and ttag is that of the other B (i?tag)- The Ai 
distribution for Bgig decaying to a CP eigenstate 



/(Ai, Qtag) 



^-\/\t\/TB 

Atb 



{1 - QtagAw + QtagM(l - 2w) 

+ (<3tag(l - 2w) + ^(1 - gtagAcj)) [5sin(AmBAi) - C cos(AmBAt)]} 



(61) 



is convolved with a resolution function TZ. The param- 
eter Qtag — +1(— 1) when the tagging meson B^^^ is a 
B'^{B'^), w is the average mistag probability, and Aw 
and /i describe the difference in mistag probability and 
the tagging efficiency asymmetry between B'^ and 5° 
mesons, respectively. The time distribution of combi- 



natorial background is assumed to have zero lifetime. 

The At resolution function TZ is the sum of three Gaus- 
sian functions representing the core, tail, and outer part 
of the distribution, weighted by the At error for the core 
and the tail: 
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TZ{At, (JAt) — fcoicG{At, I'corcO'At, CTcorcfAi) + /tailG(Ai, f tail^At , CTtailCTAt ) + foutG{At, 

^OUt ; "^OUt ) 

I 



(62) 



where G(At, cta*; i^, c) is a Gaussian distribution with 
mean and standard deviation cr, and / is the cor- 
responding fraction. We have verified in simulation 
that the parameters of the resolution function for sig- 
nal events are compatible with those obtained from the 
Bflav sample, a data sample of fully-reconstructed 
jr)(*)-7r+/p+/a+ decays. Therefore we use the i?flav pa- 
rameters for better precision. The background At dis- 
tribution is parameterized by the GP-asymmetric PDF 
/(At, ctag) = 1 ± .4bkgd(ctag), convolvcd with the reso- 
lution function. The parameters of the background At 
PDF are determined in the fit to data from sidebands in 

TOES- 

The detailed description of the treatment of the S and 
C terms of different contributing amplitudes is given in 
Section [il Dl Eq. is applicable to the time evolution 
of each of the five components in the angular distribution, 
three longitudinal {Sjo and Cjo for J = 0, 1, and 2) and 
two transverse {Sjt and Cjt for J — 1 and 2). All five 
Sjo and Sjt parameters are expressed in terms of AipQo 
and the other polarization and CP parameters entering 

the B ^ ^ ipK^TT^ PDF description, while the five Cjo 
and Cjt parameters are expressed through other polar- 
ization and CP parameters only, as shown in Eqs. ([461- 



For the combinatorial background we establish the 
functional forms and initial parameter values of the PDFs 
with data from sidebands in tties or AE. We then 
refine the main background parameters (excluding the 
resonance-mass central values and widths) by allowing 
them to vary in the final fit so that they are determined 
by the full data sample. Overall, there are 51 free back- 
ground parameters in the joint fit. 



C. Analysis validation 

We validate the analysis selection and fit performance 
with a number of cross-check analyses. To test the treat- 
ment of combinatorial background in the PDF, we per- 
form fits on the data collected below the T{4S) reso- 
nance, on GEANT-based [S^ MC simulation of about 
three times the statistics of the data sample for both qq 
production (continuum) and generic T(4S') — > BB de- 
cays. We also test the contribution of several dozen ex- 
clusive B meson decays which could potentially mimic 
the signal with statistics of more than an order of magni- 
tude greater than their expectation. No significant bias 
in the background treatment was found. Systematic un- 
certainties associated with this treatment are discussed 
in the next subsection. 

To test the signal PDF parameterization and the over- 



all fit performance, we generate a large number of MC 
experiments, each one representing a statistically inde- 
pendent modeling of the fit to the data. The signal 
events are taken from the generated MC samples, while 
background is generated from the PDF with the total 
sample size corresponding to the on-resonance data sam- 
ple. We embed signal-like events according to expecta- 
tion. We find the results of the MC experiments to be in 
good agreement with the expectations and the error es- 
timates to be correct. In Fig. [7|we show examples of the 



fitted 



pgoncratod-j^^l'^,^ distributions, where x denotes 



one of the signal parameters. The means and widths of 
all these distributions are within about 5% of the ex- 
pected values of zero and one, which results in negligibly 
small uncertainty in the fit result. 




v^rpi ^rpi 



FIG. 7: Distributions of (a;^"'^'^ - xS<="°'^="<='^)/cr(x) for a large 
number of generated MC experiments, wiiere a Gaussian fit is 
superimposed, and x denotes signal parameters; (a) /li, (b) 
(jixi, (c) A°cpi, and (d) A(^||i. 



At all analysis development steps we used a "blind" 
technique. Numerical values of all physics parameters 
Tisigj, Acpj, A(f>oQ, and C were kept hidden until the 
analysis method and tools were decided, validated, and 
fixed. Consistency between the blind fit results and pre- 
diction from generated samples for the likelihood C and 
error values were used to judge the goodness of fit. We 
also examined the projection of data and likelihood fit 
PDFs onto individual variables with the enhancement of 
either signal or background to judge consistency. De- 
tailed study of the systematic errors did not reveal any 
uncertainties which would exceed statistical errors. 
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D. Systematic uncertainties 

In Tables IVH IVIH and IVIIII we summarize the domi- 
nant sources of systematic errors in our measurements. 
In the measurement of the branching fractions, we tab- 
ulate separately the multiplicative errors on selection ef- 
ficiency in Table IVII Measurement of all parameters 
suffers from uncertainties in the fit model which are 
discussed below and in Table IVIII One additional er- 
ror in the branching fraction measurement is the error 
on the number of B'^B^ mesons produced and is esti- 
mated to be 1.1%, where we assume equal decay rates of 
T{AS) B°B° and B+B'. While most of the errors are 

dominated by uncertainties in the i3 ° — > ipK^jr^ de- 
cay mode, the measurement of AipQo has additional sys- 
tematic errors unique to the B^ — > LpKgn^ decay mode. 
Therefore all systematic errors on A0oo are quoted in 
Table Ivml 

The systematic errors in the efficiency are typically due 
to imperfect MC simulation and they are obtained from 
independent studies, such as control samples. They affect 
the errors in the branching fraction, but do not change 
the significance of the signal yield. From a study of abso- 
lute tracking efficiency, we evaluate the corrections to the 
track finding efficiencies, resulting in a systematic error of 
0.5% per track and the total error of 2.0% for four tracks. 
The error due to particle identification is about 2% and 
is dominated by the kaon selection requirements. Par- 
ticle identification performance has been validated with 
high-statistics data and MC control samples, such as D*- 
tagged D Ktt decays. The Kg selection efficiency sys- 
tematic uncertainty is taken from an inclusive Kg control 
sample study, giving a total error of 3.5%. The tt" recon- 
struction efficiency error is estimated to be 3.0% from a 
study of r decays to modes with tt" mesons. 

The reconstruction efficiency has a weak dependence 
on the fraction of longitudinal polarization due to a 
non-uniform acceptance function of the helicity angles. 
Therefore, we use the measured value of the polarization 
when computing the efficiency. The uncertainty in this 
measurement translates into a systematic error in the 
branching fraction. Several requirements on the multi- 
hadronic final state, minimum number of charged tracks, 
event-shape, and vertex requirements result in a few per- 
cent errors. Other errors come from the uncertainty in 
the daughter branching fractions for ip —^ K^K~ and 
K* Ktt. All these errors are summarized in Table IVTl 

When we perform the ML fit we make certain as- 
sumptions about the signal and background distribu- 
tions. Most background parameters are allowed to vary 
in the fit, but we constrain most B decay parameters to 
the expectations based on MC and control samples. In 
order to account for the resulting uncertainty, we vary 
the parameters within their errors, taking into account 
correlations among the parameters. We obtain ttieSj ^E, 
and uncertainties from the control samples discussed 
in Section IIIII The invariant mass uncertainties incor- 



porate errors on the resonance parameters as quoted in 
Table |IV] and Ref. [ll| . We take into account resolution 
in the Ktt and K~^K~ invariant masses with the corre- 
sponding errors on the absolute values of 1 MeV and 0.3 
MeV, respectively. 

We separate a special class of PDF uncertainties for the 
helicity angles, due to the acceptance function. In addi- 
tion to statistical errors in the MC sample, we consider 
the momentum-dependent uncertainty of the tracking ef- 
ficiency. The main effect is on the curvature of the accep- 
tance function shown in Fig. [5] due to a strong correlation 
between the momentum of a track and the value of the 
helicity angle. Moreover, in order to study the effects of 
charge asymmetry in angular distributions, we apply the 
acceptance correction independently to only B or only B 
decay subsamples. The largest deviation is taken as the 
"acceptance" systematic error quoted in Table IVlIl 

To estimate the effect of the B meson decays which 
could mimic signal, we study a full GEANT4-based MC 
simulation of the T(4S') BB events. We embed the 
categories of events which may have AE and jties dis- 
tributions similar to signal into the data sample and ob- 
serve variations of the fit results which we take as sys- 
tematic errors. The nonresonant contribution is taken 
into account naturally in the fit with both KK and 
Ktt contributions allowed to vary. The former is mod- 
eled as _B" foK*^ and the latter is a part of the 
5-wave Ktt parameterization. Interference effects are 
studied separately. We also take into account the un- 
certainty in the shape of the K^ K~ invariant mass dis- 
tribution. The default parameterization assumes the 
B° foK*{892)° decay and we vary it to the phase- 
space B° K+K-K*{892)° distribution. We constrain 
the number of (pK* (892)'^ events contributing to 

the higher Ktt invariant mass range based on the mea- 
sured branching fraction, but we also vary this number 
according to the branching fraction errors. These errors 
are quoted as "B-bkgd" in Table EIIl 

The selected signal B'^ — ipK*^ events contain a small 
fraction of incorrectly reconstructed candidates. Misre- 
construction occurs when at least one candidate track 
belongs to the decay products of the other B from the 
T{AS) decay, which happens in about 5% of the cases in 

the K *^ ^ K^TT^ decay. The distributions that show 
peaks for correctly reconstructed events have substantial 
tails, with large uncertainties in MC simulation, when 
misreconstructed events are included. These tails and 
incorrect angular dependence would reduce the power 
of the distributions to discriminate between the back- 
ground and the collection of correctly and incorrectly re- 
constructed events. We choose, therefore, to represent 
only the correctly reconstructed candidates in the signal 
PDF, and to calculate the reconstruction efficiency with 
both the correctly reconstructed and misreconstructed 
MC events. Fitting the generated samples to determine 
the number of correctly reconstructed candidates has an 
efficiency close to 100% even though a few percent of 
selected candidates are identified as background. We 
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TABLE VI: Systematic errors (%) in reconstruction efficiency evaluation. The total errors combine the two subchannels 

according to their weight and are dominated by the K^n^ channel. We separate the ip{K'K)llP ^'^^ <^ifo (1430)" modes to 
account for different errors in the daughter branching fractions. See text for details. 





ipK*{m2f 


(^/s:2*(1430)° 


ipyi^-n ) 
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track nndmg 


2.0 


1.0 


2.0 


1.0 


2.0 


1.0 


2.0 




1.0 




2.1 


1.1 


2.1 


1.1 


2.1 


1.1 


2.1 




1.1 


Kg selection 




3.5 




3.5 




3.5 






3.5 


Tr" selection 




3.0 




3.0 




O.U 






3.0 


MC statistics 


0.2 


0.3 


0.2 


0.3 


0.3 


0.4 


0.3 




0.4 


polarization 


0.2 


1.1 


0.2 


1.7 












event selection 


1.0 




1.0 




1.0 






1.0 




thrust angle 6t 


1.0 




1.0 




1.0 






1.0 




vertex requirement 


2.0 




2.0 




2.0 






2.0 




branching fraction 


1.2 




1.2 




1.2 






1.2 




K* branching fraction 


0.0 




2.4 




0.0 






10.8 




K% branching fraction 




0.1 




0.1 




0.1 






0.1 



total 4.0 4.7 4.0 11.5 



TABLE VII: Systematic errors in the measurement of the three signal yields (%) and other signal parameters (absolute values), 
excluding the l^(j)m measurement. Uncertainties due to parameterization, acceptance function modeling, B-background (B- 
bkgd), fit response, interference of the (KK) final states (intcrf.), charge asymmetry in reconstruction, assumptions about the 
unconstrained CP asynunctries Acp2 a-^d A4>±2 {CP asym.), and the total errors are quoted. The errors are not quoted if they 
are either small or not relevant for a particular measurement. See text for details. 





PDF 
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B-bkgd 
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interf. 


charge 


CP asym. 


total 


ifiiKnYo" yield (%) 


7.1 




3.4 


2.7 








8.3 


^K*{892f yield (%) 


2.3 




1.9 


1.3 


1.8 






3.7 


cpKi{U30)° yield (%) 


3.3 




0.8 


2.8 








4.4 


Acpo 


0.036 


0.002 


0.048 


0.001 




0.020 


0.008 


0.064 


hi 


0.002 


0.002 


0.005 


0.007 


0.010 






0.013 


f±i 


0.001 


0.001 


0.006 


0.004 


0.010 






0.013 


hi 


0.007 


0.001 


0.010 


0.017 


0.078 






0.081 


(l>±i 


0.005 


0.001 


0.010 


0.010 


0.084 






0.085 


Acpi 


0.007 


0.001 


0.008 


0.012 


0.018 


0.020 




0.031 


■^CPl 


0.003 


0.003 


0.013 


0.005 


0.019 






0.024 


Acpi 


0.007 


0.004 


0.022 


0.024 


0.052 






0.062 


A0II1 


0.008 


0.001 


0.009 


0.010 
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0.080 


A0X1 


0.007 


0.001 


0.007 


0.016 


0.081 






0.083 


Soi 


0.030 
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0.005 


0.003 


0.080 






0.086 


ASoi 


0.008 


0.001 


0.006 


0.006 


0.080 






0.081 


fL2 


0.006 
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0.016 
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0.004 


0.037 


f±2 


0.001 


0.001 


0.004 


0.031 






0.003 


0.031 


0||2 


0.051 


0.002 


0.021 


0.029 






0.012 


0.063 


AcP2 


0.037 


0.002 


0.029 


0.001 




0.020 


0.005 


0.051 


•f^CP2 


0.012 


0.001 


0.005 


0.002 






0.003 


0.014 


A0II2 


0.027 


0.002 


0.088 


0.012 






0.011 


0.093 


502 


0.117 


0.004 


0.062 


0.009 






0.006 


0.133 


A5o2 


0.012 


0.004 


0.057 


0.009 






0.007 


0.059 
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TABLE VIII: Systematic errors (absolute values) in the mea- 
surement of A(j)oo. See Table IVIll and text for details. 



TB, Ams 0.001 

sin2/3 measurement 0.015 

signal At resolution 0.016 

mistag differences 0.019 

2 scale -|- boost 0.002 

beam spot 0.010 

SVT alignment 0.001 

tag-side interference 0.002 
background resolution and asymmetry 0.006 

B-background 0.024 

PDF 0.009 

acceptance 0.004 

fit 0.002 

CP-asym. 0.003 

total 0.041 



account for this with a systematic error taken as half 
the fraction of candidates identified as background and 
quoted as the "fit" entry in Table IVIII Similarly, we 
obtain systematic errors on other parameters from the 
largest deviation from the expectation. This includes a 
potential bias from the finite resolution in the helicity 
angle measurement and a possible dilution due to the 
presence of the misreconstructcd component. 

As we discuss below, a substantial 5° foKn contri- 
bution is found in the lower K-k mass range, correspond- 
ing to either /o-ftr*(892)° decays, or any other con- 
tribution with a broad K~^K~ invariant mass distribu- 
tion, either resonant or nonresonant. The uncertainties 
due to rrij^j^ interference are estimated with the samples 
generated according to the observed K^K^ intensity and 
with various interference phases analogous to i5oj in K-k. 
These are the dominant systematic errors for the C pa- 
rameters of the B^ — > LpK*{S92Y' decay. No significant 
B^ foKn contribution is observed in the higher Ktt 
mass range. 

The charge bias uncertainty affects only the relative 
yields of B and B events. We assign a systematic error 
of 2%, which accounts mostly for a possible asymm etry 
in the reconstruction of a charged kaon from a K*^ [33| ■ 
Overall charge asymmetry has a negligible effect on the 
angular asymmetry parameters, while the angular depen- 
dence of the charge asymmetry is tested with the flavor- 
dependent acceptance function discussed above. 

There are still two CP parameters, Aqp2 and A(f)±2, 

that are not measured in the ip'^K^llASO)^ 
(fK^TT^ decay. We assume zero asymmetry as the most 
likely value and vary them within ±0.2 for the direct- 
CP asymmetry and ±0.5 rad for the phase asymmetry. 
All the errors on the fit parameters are summarized in 
Table EUl 



For the time-dependent measurement of Ac^qq, the 
variations are quoted in Table IVTIIl We vary the B'^ life- 
time and Am_B within their errors We include the 
error of the (3 measurement from Eq. We use the 

results of the sin2/3 analysis ^3j] to estimate the system- 
atic errors related to signal parameters, such as At signal 
resolution and mistag differences, detector effects (z scale 
and boost, beamspot, SVT alignment uncertainties), and 
the tag-side interference. Background CP-asymmetry 
and resolution parameters are determined from the side- 
band data and then constrained in the fit. The con- 
strained parameters are varied according to their errors 
and added in quadrature to compute the systematic er- 
rors. 



V. RESULTS 

We observe a non-zero yield with more than 10 a sig- 
nificance, including systematic uncertainties, in each of 
the three i3° — > (pK*° decay modes. In Figs. IMTUl we 
show projections onto the discriminating variables. For 
illustration, the signal fraction is enhanced with a re- 
quirement on the signal-to-background probability ratio 
to be greater than a value within the range (0.85 - 0.95), 
calculated with the plotted variable excluded. This re- 
quirement is at least 50% efficient for the signal events. 

In Tables ITXl and 1x1 the rigigj, Acpj, A^oo, and C = 

{/l,7, f±J, (l)\\J, (f>±.J, Sqj, ^cPj, Acpj, A(^||j, A(/)_lj, 

ASqj} parameters of the B° ipK*{892)°, ipK^iUSOf, 
and (p{K'!r)Q'^ decays are shown. The three quantities 
0±2, Aqp2' ^^'^ A(j>±2, which characterize parity-odd 
transverse amplitude in the vector-tensor decay, are not 
measured because f±2 is found to be consistent with zero. 

The computed significance of the yield is more than 
24 cr for 5° tpK* {892)° and 11 cr for B° ip{K7r)^°. 
Given the convincing presence of the 5- wave {Ktt)q'^ con- 
tribution, we rely on the interference terms to resolve the 
phase ambiguities. In the lower mK-rr range the yield of 
the (p{Ktt)q^ contribution is 751^° events with a statis- 
tical significance of 9 cr, including the interference term. 
From the measurements of the higher mK-ir range we cal- 
culate the contribution of (^(iir7r)5° to the lower mass 
range to be about 61 events. This is consistent with the 
above result within 1 cr. The dependence of the inter- 
ference on the Ktt invariant mass [2^, |23| allows us to 
reject the other solution near (27r — tt — (j)±i) relative 
to that in Table |X] for the B" ipK* {892f decay with 
a statistical significance of 6.5 cr (which becomes 5.4 cr 
when systematic uncertainties are included). We also re- 
solve this ambiguity with a statistical significance of more 
than 4 cr with the i?° or i?° decays independently. Fig. [TT] 
shows the = — 21n(£/£niax) scan plots for 0|| and (l>±, 
where we illustrate how the phase ambiguity is resolved. 
For comparison, we show the result of the fit where in- 
terference is not taken into account and no sensitivity to 
resolve the ambiguity is present. The significance of the 
deviations of and from tt is 5.4 cr (4.5 cr) and 6.1 cr 
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TABLE IX: Analysis results: the reconstruction efBciency Ercco obtained from MC simulation; the total efficiency e, including 
the daughter branching fractions Ts^; the number of signal events Wsigj; statistical significance (<S) of the signal; the branching 
fraction Bj\ and the flavor asymmetry Acpj- The branching fraction B{B^ ip{KTv)o^) refers to the coherent sum \ArcB + 
j4non-resP of resouaut and nonresonant = 0''' Ktt components and is quoted for mxir < 1-6 GeV, while the B{B'^ — » 
(pi("o (1430)") is derived from it by integrating separately the B-W formula of the resonant |ylrcs|^ Kit component without rriKir 
restriction. The systematic errors are quoted last. 



mode 


Ercco (%) 




Jisigj (events) 


S (a) 


Bj (10-'^) 


Acpj 












3.9 ±0.5 ±0.6 






23.2 ±0.9 
11.7 ±0.4 


8.3 ±0.3 
7.6 ±0.3 
0.66 ±0.03 


172 ± 24 ± 14 
158 ± 22 ± 13 
14±2± 1 


11 


4.3 ±0.6 ±0.4 


±0.20 ±0.14 ±0.06 


V5-ft:*(892)° 


33.7 ± 1.3 

13.8 ±0.5 


11.9 ±0.4 
11.1 ±0.4 
0.78 ±0.03 


535 ± 28 ± 20 
500 ± 26 ± 19 
35 ± 2 ± 1 
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9.7 ±0.5 ±0.5 


±0.01 ±0.06 ±0.03 


957^1(1430)° 


26.7 ± 1.0 
8.7 ±0.3 


4.7 ±0.2 
4.4 ±0.2 
0.25 ±0.01 


167 ±21 ±8 
158 ± 20 ± 7 
9± 1 ± 1 


11 


7.5 ±0.9 ±0.5 


-0.08 ±0.12 ±0.05 



(S.Oct), respectively (including systcmatics in parenthe- 
ses). 

Projections onto Tii and Ti.2 in the lower rriKTr range 
in Figs. [H (a) and (b) show sizable contributions of 
both cos^ 9 (longitudinal) and sin^ 9 (transverse) com- 
ponents. These two plots emphasize /n {Hi, 7^2, ^) 
and /21 {Hi, H2, $) angular terms in Eq. ([5]). Similarly, 
projections onto Hi and H2 in the higher mpcir range 
in Fig. [TU] show predominant longitudinal polarization 
with sizable /12 {Hi, H2, ^) and fio {Hi, H2, <&) contri- 
butions, while the /22 {Hi, H2, $) contribution is small. 

In order to illustrate /31 {Hi, H2, $) and 
/41 {Hi, H2, $) angular distributions, we project 
onto the angle — $ for the B decays and <& for B decays. 
This procedure takes into account the change of sign for 
the odd components with P-wave amplitude {Ai±). This 
$ angle projection is sensitive only to the constant term 
and to cos(2$) and sin(2<i>) terms. The fact that the 
double sine and cosine contributions are small in Fig. [9] 
(c) tells us that both (|Ai|| p - |Ai_Lp) and '<sm{Ai±Al^^) 
are relatively small, in agreement with the fit results. 
That is, the values of (1 — Jli — 2/_li) and (^p — (j)±i) 
are small; see a^j and a^j in Eqs. pT|) and (dH). 

In order to emphasize /51 {Hi, H2, $) and 
f6i{Hi,H2,^) angular terms in Eq. or cosi* 

and sin$ distributions, we show the difference be- 
tween the above $ angle projections for events with 
Hi X H2 > and with Hi x H2 < 0. This gives us 
contributions to a^j and a^j in Eqs. (|13p and (|14p . 
while background and all other signal contributions 
cancel. This projection is shown in Fig. [9] (d), where we 
see good agreement between the data and the fit results. 
This plot indicates a sizable cos $ component due to 
5Re(j4i|| AJq) and an asymmetric sin<I> component due 



to '^Tn{Ai±AlQ). The latter asymmetry is visible and 
reflects the presence of a strong phase with a statistical 
significance of 6.4 ct. 

In order to emphasize CP asymmetries in the an- 
gular distributions we make two other similar projec- 
tions onto the $ angle, but now we plot the differ- 
ence between the B and B decays. Figs. |9] (e) and 
(f) show distributions of (Ns>— Ng^+NB<— N^^) and 
(N_B>— Ng^— Ns<+Ng^), where Nb> denotes the num- 
ber of B decays with H1H2 > 0, N|j^ is the number 
of B decays with H1H2 > 0, Nb< is the number of B 
decays with H1H2 < 0, and N^^ is the number of B 
decays with H1H2 < 0. In all cases we project onto — $ 
for the B decays and onto <i> for the B decays. Fig. [9] 
(e) is sensitive to CP asymmetries in the cos(2$) and 
sin(2$) terms in Eqs. ^ and dH]), while Fig. [9] (f) is 
sensitive to CP asymmetries in the cos $ and sin $ terms 
in Eqs. (fT3|) and (fT4|) . In particular, there is a hint of a 
sine wave contribution in Fig. [5](f), corresponding to the 
non-zero measurement of A^j^i — +0.21 ± 0.13 ± 0.08, 
though not significant enough to constitute evidence for 
CP violation. 

In Table |X] we summarize the correlation among the 
primary fit parameters. There are two large correla- 
tion effects of approximately 50% evident for the i?" 
(pK* {892)^ decay in the fit, that is between the longitu- 
dinal and transverse fractions {Jli and f±i) and between 
the two phases (0|[i and 4i±i). In Fig.[T2]we show the like- 
lihood function contour plots for the above pairs of corre- 
lated observables as well as for and f±2- Fig.fOlshows 
the = — 21n(£/£max) distributions for the CP viola- 
tion phase parameters ASxj and A(f>xj, where x stands 
for either ±, ||, or 0. 

The — > JqK*^ category accounts for final states 
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TABLE X: Summary of the results, see Table |II] for definition of the parameters. The branching fractions Bj and flavor 
asymmetries Acpj are quoted from Table HXl The systematic errors are quoted last. The dominant fit correlation coefficients 
(p) are presented for the (pK* (892)'^ and (^7^2(1430)" modes where we show correlations of So with (f)\\/(f)± and of ASo with 
A0||/A(;!>x. 



parameter 


J = 


J = 1 


P 


V9_K"2(1430)° 
J = 2 


P 


Bj (10-'^) 


3.9 ±0.5 ±0.6 


9.7 ±0.5 ±0.5 




7.5 ±0.9 ±0.5 




fLJ 

.f±J 




0.494 ± 0.034 ± 0.013 
0.212 ±0.032 ±0.013 


}-48% 


0-90llo;o58 ± 0.037 
0.002to ± 0.031 


}-15% 


<t>\\J (■'^'i) 
(j)^j (rad) 




2.40 ±0.13 ±0.08 
2.35 ±0.13 ±0.09 


} 62% 


3.96 ±0.38 ±0.06 


- 


ooj (Taa) 




9 SO -1-0 1 -1- n HQ 


o4 /o / zo /o 


o/ji _|_niQ_|_nio 
o.4i zt U.lo zt U.lo 


iy /o 


Acpj 


+0.20 ± 0.14 ± 0.06 


±0.01 ± 0.06 ± 0.03 




—0.08 ± 0.12 ± 0.05 




Acpj 
Acpj 




±0.01 ±0.07 ±0.02 
-0.04 ±0.15 ±0.06 


}-47% 


-0.05 ±0.06 ±0.01 




A(^l|j (rad) 
A(f>±j (rad) 




±0.22 ±0.12 ±0.08 
±0.21 ±0.13 ±0.08 


} 62% 


-1.00 ±0.38 ±0.09 




ASoj (rad) 




±0.27 ±0.14 ±0.08 


35%/24% 


±0.11 ±0.13 ±0.06 


16% 


A<^oo (rad) 


0.28 ±0.42 ±0.04 











with K^K^ from either /o, ao, or any other broad 
K~^K~ contribution under the ip. Its yield is consistent 
with zero in the higher rriKTr range and is 84± 19 events in 
the lower mxTr range. Due to uncertainties in the nature 
of this contribution, we do not calculate its branching 
fraction, but include it in the evaluation of systematic 
uncertainties in other parameters as discussed above. 

In Fig.[T3]wc illustrate the effect of interference in the 
Kit invariant mass. In the vector-scalar Kn interference 
(lower mass range) the interference term is linear in Tii, 
creating a forward-backward asymmetry. However, due 
to variation of the B-W phase this effect cancels when 
integrated over the mxTr range. For i5o — tt we expect 
the coefficient in front of Tii to be positive above rriKir — 
0.896 GeV and negative below this value. Thus, we create 
Fig. [13] (a) to emphasize this effect. If Sqi ~ 0, the sign 
of the forward-backward asymmetry would be reversed 
and we would see more events on the left as opposed to 
the right. Thus Sqi ~ tt is preferred. We note that this 
plot has only partial information about the interference 
while the multi-dimensional fit uses the full information 
to extract the result. 

Figure [13] (b) shows a similar effect in the tensor-scalar 
Kit interference (higher mass range). In this case inter- 
ference could either enhance events in the middle of the 
Til distribution and deplete them at the edges, or the 
other way around. Fig. [TJ] (b) indeed shows significant 



improvement in the Tii parameterization with the inclu- 
sion of interference. It also corresponds to the observed 
value Sq2 — tt. 

Because of the low significance of our measurements 
of /||2 = (1 - /l2 - /±2) (1.9a) and f^2 (Oa) in the 
B'^ (1430)° decay we have insufficient information 

to constrain 0||2 at higher significance or to measure (f>±2, 
Aqp2, or A(j>±2, which we constrain to zero in the fit. 

Finally, we fit a single parameter A0oo in the time 
evolution after combining all available i?° — ^ (pKgTr'^ 
charmless final states, which are dominated by spin-0, 
1, and 2 (Ktt) combinations. The distribution of the 
time difference At and the time-dependent asymmetry 
are shown in Fig. 1151 The parameter Acjjoa is measured 
to be 0.28 ± 0.42 ± 0.04, as shown in Table |Xj 



VI. CONCLUSION 

In conclusion, we have performed an amplitude anal- 
ysis and searched for CP violation in the angular dis- 
tribution of i?" ipK*^ decays with tensor, vector, 
and scalar K*'^ mesons. Our results are summarized 
in Tables IIXI and 1X1 a nd supersede corresponding mea- 
surements in Ref. [20|. In this analysis we employ sev- 
eral novel techniques for CP-violation and polarization 
measurements in the study of a single _B-decay topol- 
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FIG. 8: Projections onto the variables mES, Ai?, raK-rr, 
and m^j^ for the signal B ipK^-K^ (left) and B — > tfiB^n'^ 
(right) candidates. Data distributions are shown with a re- 
quirement on the signal-to-background probability ratio cal- 
culated with the plotted variable excluded. The solid (dashed) 
hues show the signal-plus-background (background) PDF pro- 
jections. 



FIG. 9: Projections onto the variable Tii for the lower tuktv 
range in (a) and Ti.2 in (b). Projections onto Q x "I> for the 
lower rriK-K range in (c)-(f) where Q changes sign between the 
B decays and the B decays. The distributions are shown 
for the signal —^ (pK* (892)" candidates following the solid 
(dashed) line definitions in Fig. [8] The $ angle projections are 
shown for different combinations of event yields with certain 
requirements on the B flavor and Hi x Ti.2 product signs, as 
discussed in the text. The D^j-meson veto causes the sharp 
acceptance dips near TCi = 0.8 seen in (a). 





FIG. 10: Same as Fig. |9] (a,b), but for the signal B° 
(^ifj (1430)" and lp{K-k)o^ candidates combined. The D^^y 
meson veto causes the sharp acceptance dips near Tii — 0.8 
seen in (a). 



ogy (p{Ktt). We use the time-evolution of the 

(pKgw'^ channel to extract the (7P-violating phase 
difference A^qq = 0.28 ± 0.42 ± 0.04 between the B and 
B decay amplitudes. We use the dependence on the Ktt 
invariant mass of the interference between the scalar and 
vector, or scalar and tensor components to resolve dis- 
crete ambiguities of both the strong and weak phases. 
Twelve parameters are measured for the vector-vector 
decay, nine parameters for the vector-tensor decay, and 
three parameters for the vector-scalar decay, including 



the branching fractions, CP-violation parameters, and 
parameters sensitive to final-state interactions. 

The {V-A) structure of the weak interaction and the s- 
quark spin-flip suppression in the process shown in Fig. [1] 
suggest l^ljol > \Aj+i \ > |yl.7_i| [ll,[3. The relatively 
smaU value of /li = 0.494 ± 0.034 ± 0.013 and the rela- 
tively large value of f±i ^ 0.212 ± 0.032 ± 0.013 in the 
vector-vector decay remain a puzzle. The naive expecta- 
tion is that the {V-A) nature of the weak decays requires 
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FIG. 11: Scan of X = -2 ln(£//:max) as a function of ^yi and 
0X1 for the (pK^n^ decays in the lower itiktv range, where the 
filled circles are the results with the interference term, and 
the open circles without the interference term. Two discrete 
solutions are visible in the case without interference, while 
this ambiguity is resolved with the interference term. The 
values of A^yi and A<j!>xi have been constrained in the range 
(—0.5, 0.5) in order to reject ambiguities with larger values of 
A(/!>||i and A^xi. 



that an anti-quark originating from the b — s- qW~^ de- 
cay be produced in helicity state This argument 
applies to the penguin loop which is a purely weak tran- 
sition, b ^ s, with the double-M^ coupling in the Stan- 
dard Model. The s anti-quark can couple to the s quark 
(see Fig. [Ij to produce the (p state with helicity of either 
A = or A = +1, but not A = —1. However, the K* 
state should have the same helicity as the cp due to an- 
gular momentum conservation. The A = +1 state is not 
allowed in this case because both s and s quarks would 
have helicity in violation of helicity conservation in 
the vector coupling, g ss. 

The spin-flip can alter both of the above requirements, 
but its suppression factor is of order ~ my /ttib for each 
flip, where my is the mass of the ip or A'*" mesons. Thus, 
we arrive at the expectation \Ajo\ ^ |^j+i| ^ 
or |Ajo| > \Aj±\ and Aj± ~ Aj\\, where Aj+i is sup- 
pressed by one spin flip, while Aj^i is suppressed by 
two spin flips. New physics could have different inter- 
actions, alter the spin-helicity expectations and result 
in a large fraction of transverse polarization. Alterna- 
tively, strong interaction effects might change this expec- 



+0.046 
-0.058 



± 0.037 in 



tation [161. The value of /l2 = 0.901 
vector-tensor decays is not compatible with that mea- 
sured in vector-vector decays, while it is compatible with 
the expectation from the spin-flip analysis above. This 
points to a unique role for the spin-1 particle recoiling 
against the ip in the B ipK* polarization puzzle. 

In the ipK* {892)° decay we obtain the solu- 

tion ~ 0_Li without discrete ambiguities. Com- 
bined with the approximate solution f^i ~ 1/2 and 
f±i ~ (1 — /ii)/2, this results in the approximate de- 
cay amplitude hierarchy |Aio| ~ |^i+i| ^ l^i-i| (and 
\Aio\ ~ \Ai-i\ > |^i+i|). We flnd more than 5a (4ct) 
deviation, including systematic uncertainties, of 4)±{(f>\\) 
from either tt or zero in the ipK* (892)'^ decay, in- 

dicating the presence of flnal-state interactions (FSI) not 
accounted for in the naive factorization. The effect of FSI 
is evident in the phase shift of the cosine distribution in 



Fig. [91(d). 

From the definition in Eq. (|46p and the measurement 
of A^ooj we determine the parameter 



sin(2/3eff ) = sin(2^ + 2A0oo) = 0.97: 



-0.03 
-0.52 I 



(63) 



as measured with the B'^ ipKq{1A30)^ decay. Our 
measurements of eleven CP- violation parameters rule out 
a significant part of the physical region and are consistent 
with no CP violation in the direct decay, but are consis- 
tent with the measurement of sin(2/3off) in Eq. ([63)) . The 
current precision on sin(2/3off) is still statistics-limited, 
although we constrain sin(2/3eff) > 0.15 at the 90% confi- 
dence level, which is consistent with the Standard Model 
CP violation due to B^-B^ mixing. This analysis pro- 
vides techniques for future experiments to extract this 
parameter from the B^ ipKn decays. 

Other significant non-zero CP-violation parameters 
would indicate the presence of new amplitudes with dif- 
ferent weak phases. The parameters and A(/)||j 
are particularly interesting due to their sensitivity to the 
weak phases of the amplitudes without hadronic uncer- 
tainties [13], such as the relative weak phases of Aj+i 
and A JO, while the CP- violation parameter ASqj repre- 
sents potential differences of weak phases among decay 
modes. 

We note that the measurement of sin(2/3eff ) in Eq. ([63[) 
is not the primary result of this analysis, but only an 
interpretation of the A0oo measurement. Equivalently, 
there could be six other effective sin(2/3) measurements 
as shown in Eqs. ([T7H15[) . However, all of them would be 
highly correlated due to the same dominant uncertainty 
coming from A(f)oo- Rather than give them all here, we 
provide an illustration of our measurements with the fol- 
lowing differences using the results in Table 1X1 as input: 



sin(2/3-2A(5oi)-sin(2/3) 



_n 40-1-0. 26 

"•^^-0.34 ; 



(64) 



sin(2/3 - 2A^||i) - sin(2/3) = -0.32l° , (65) 

sin(2/3 - 2A0ii) - sin(2/3) = -0.30l^;3^ , (66) 
sin(2/3 - 2A0^i) - sin(2/3 - 2A0||i) = 0.02 ± 0.23 , (67) 
sin(2/3 - 2A(5o2) - sin(2/3) = -0.101E;;29 . (68) 

Systematic uncertainties are included in the errors 
quoted in Eqs. ([63H68| . 

Taking the example in Eq. ([S71) . we see that because 
of the positive correlation between A(f>±i and A(/)||]^, we 
achieve a precision of ±0.23 on the measurement of the 
difference between values of sin(2/9off) from the parity- 
odd {Ai±) and parity-even (^i||) decay amplitudes. This 
precision is significantly better than that of the measure- 
ment of sin(2/3off) itself because of the cancellation of 
common uncertainties. A significant deviation from zero 
would indicate a CP-violating contribution to either the 
parity-odd or parity-even amplitude but not the other. 
A similar comparison would be the values of sin(2/?off) 
measured in B ^ rj'K and B ipK decays. This mea- 
surement of the sin(2/3cff) difference with the parity-odd 
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FIG. 12: Contours corresponding to the unit intervals ol \/~2A \nC for polarization f±j and /lj (a) and phase (f)±i and 
(b) measurements. Diagonal dashed lines f±j = (1 — fLj)/2 and — correspond to |j4j+i| ^ In (a), the solid 

(dashed) contours are the results for J = 1 (J = 2). In (b) the (7r,7r) point is indicated by the crossed dashed lines. 




-0.2 0.2 0.4 0.6 0.8 
A5o, (rad) 



.S 10 



(C) 



-0.2 0.2 0.4 0.6 
A(^„, (rad) 




-0.2 0.2 0.4 0.6 
A(^, (rad) 




-0.2 0.2 0.4 0.6 0.8 
A5o2 (rad) 



t 2 



(d) 



-1.5 



-0.5 



A(^„, (rad) 




1 
A6,(rad) 



FIG. 13: Scan of = -2\n{C/Cm^^) as a function of (a) 
A(5oi, (b) A(So2, (c) A<ji||i, (d) A^ya, (e) A^xi, and (f) A^oo- 



and parity-even amplitudes is possible with the angular 
analysis alone without any time-dependent measurement. 

Among other results in this analysis, we note the sig- 
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FIG. 14: Projections onto the variable ±7ii. (a) low tuktv 
mass range, where we project the data onto Tii for mKn > 
0.896 GeV and onto -Hi for m^^ < 0.896 GeV. The points 
with error bars show data and the histogram corresponds to 
the results of the MC generated with the observed polariza- 
tion parameters, (b) high rriKTv mass range, where we project 
the data onto Hi. Points with error bars represent the data, 
while the solid line represents the PDF projection with the 
interference term and the dashed line without the interfer- 
ence term included. The D^j-meson veto causes the sharp 



acceptance dips near Ti.i = 0.8. 



nificant yield of events (more than 5 a statistical signif- 
icance) in the category B° -> {K+K-)K*{892)°, where 
(K'^K^) reflects an S'-wave contribution, which could 
be an /o or ag meson, or any other scalar component. 
This decay is of the type scalar- vector. We have al- 
ready observed such decays with B° ~f (ysi^Q (1430)°, 
as discussed in this paper. Therefore it is plausible 
that the two decays are related by SU(3) symmetry. 
However, we do not report the branching fraction of 
B° {K+K-)K* (892)" because the exact nature of the 
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FIG. 15: The distribution of At for events in the signal 
region, for -Stag (a-) a-nd B'^a.g, (b) events with the fit re- 
sult overlaid. The solid (dashed) lines show the signal-plus- 
background (background) PDF projections. The asymmetry 
A{At) shown in (c) is defined in Eq. pi)) . 



process is not known, and a more detailed study together 
with B° -> (7r+7r-)iv:*(892)° is required. Nonetheless, 



interference between the B° {K+K-)K* {892)° and 
— > ipK*{892)^ decays provides a further path for re- 
lating strong and weak phases in the two processes, sim- 
ilar to the interference studies presented in this analysis. 
At present this interference is considered in the study of 
systematic uncertainties with good prospects for phase 
measurements from higher statistics experiments. 
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